Analysis of Generalized Grover's Quantum Search Algorithms 
Using Recursion Equations 
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The recursion equation analysis of Grover's quantum search algorithm presented by Biham et 
al. [PRA 60, 2742 (1999)] is generalized. It is applied to the large class of Grover's type algorithms 
in which the Hadamard transform is replaced by any other unitary transformation and the phase 
inversion is replaced by a rotation by an arbitrary angle. The time evolution of the amplitudes of 
the marked and unmarked states, for any initial complex amplitude distribution is expressed using 
first order linear difference equations. These equations are solved exactly. The solution provides the 
number of iterations T after which the probability of finding a marked state upon measurement is 
the highest, as well as the value of this probability, P ma x- Both T and P max are found to depend 
on the averages and variances of the initial amplitude distributions of the marked and unmarked 
states, but not on higher moments. 
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I. INTRODUCTION 

Grover's search algorithm Jl],[| provides a dramatic 
example of the potential speed-up offered by quantum 
computers. It also provides an excellent laboratory for 
the analysis and implementation of quantum algorithms 
in various hardware media. The problem addressed by 
Grover's algorithm can be viewed as trying to find a 
marked element in an unsorted database of size N. To 
solve this problem, a classical computer would need, on 
average, N/2 data base queries and N queries in the 
worst case. Using Grover's algorithm, a quantum com- 
puter can accomplish the same task using merely 0(y/~N) 
queries. The importance of Grover's result stems from 
the fact that it proves the enhanced power of quantum 
computers compared to classical ones for a whole class 
of oracle-based problems, for which the bound on the 
efficiency of classical algorithms is known. 

Grover's algorithm can be represented as searching 
a pre-image of an oracle-computable boolean function, 
which can only be computed forward, but whose in- 
verse cannot be directly computed. Such a function is 
F : D — > {0,1} where D is a set of N domain values 
(or states) and the pre-images of the value 1 are called 
the marked states. The problem is to identify one of the 
marked states, i.e., some v €E D such that F(v) = 1. 
Problems of this type are very common. One important 
example, from cryptography, is searching for the key K of 
the Data Encryption Standard (DES) |f], 

given a known 

plaintext P and its ciphertext C, where F = 1 if the 



pair of plaintext and ciphertext match [i.e., Ek{P) = C 
where Ek is the encryption function] and F = other- 
wise. Other examples are solutions of nondctcrministic 
polynomial time (NP) and NP-complete problems, which 
include virtually all the difficult computing problems in 
practice 

Let us assume, for simplicity that N = 2", where n is 
an integer. We introduce a register \x) = \x\ . . . x n ) of n 
qubits to be used in the computation. Grover's original 
quantum search algorithm consists of the following steps: 

1. Initialize the register to |0) = 1 . . . 00) , and ap- 
ply the Hadamard transform to obtain a uniform 
amplitude distribution. 

2. Repeat the following operation T times: 

A. Rotate the marked states by a phase of it radi- 

ans. 

B. Rotate all states by ir radians around the aver- 

age amplitude of all states. This is done by 
(i) Hadamard transforming every qubit; (ii) 
rotating the |0) state by a phase of ir radians, 
and (iii) again Hadamard transforming every 
qubit. 

3. Measure the resulting state. 

A large number of results followed Grover's discovery. 
These results include the proof || that the algorithm is 
as efficient as theoretically possible ||. A variety of ap- 
plications were developed, in which the algorithm is used 
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in the solution of other problems • Recently, exper- 

imental implementations of Grover's algorithm were con- 
structed using nuclear magnetic resonance (NMR) [ ^5|Jl6| 
as well as an optical device JItJ . 

Several generalizations of Grover's original algorithm 
have been developed. The case in which there are sev- 
eral marked states was studied in Ref. |jl8|| . Let k(t) 
[l(t)] denote the amplitude of the marked [unmarked ] 
states after t iterations of the algorithm. It was shown 
in that the amplitude of the marked states in- 
creases as: k(t) = sin[cj(t + 1/2)]/ y/r, while at the same 
time that of the unmarked states decreases as: lit) = 
cos[w(i + l/2)]/y/N — r, where u = 2 axcsm(\J r/N) and 
r is the number of marked states. For N ^> r the opti- 
mal time to measure and complete the calculation is after 
T = 0{<JN/r) iterations, when kit) is maximal. 

Recently, the algorithm was further generalized by al- 
lowing an arbitrary (but constant) unitary transforma- 
tion to take the place of the Hadamard transform in the 
original setting fll9|-pl| and an arbitrary phase rotation 
of the marked and the predefined states instead of the 
7r inversion j2S| . Another generalization was obtained by 
allowing for an arbitrary complex initial amplitude distri- 
bution, instead of the uniform initial amplitude distribu- 
tion obtained in step 1 above 23 24j . 

In this paper we analyze the time-evolution of the am- 
plitudes in the most general algorithm, using an arbitrary 
unitary transformation and phase rotations on an arbi- 
trary complex initial amplitude distribution. Using first 
order linear difference equations we obtain an exact so- 
lution for the time evolution of the amplitudes under the 
most general conditions. The solution provides the opti- 
mal number of iterations T after which the probability of 
finding a marked state upon measurement is the highest, 
as well as the value of this probability, P max . Both T and 
-Pmax are found to depend on the averages and variances 
of the initial amplitude distributions of the marked and 
unmarked states, but not on higher moments. 

The paper is organized as follows. In Sec. II we present 
the generalized algorithm. The analysis based on recur- 
sion equations is given in Sec. III. The results are dis- 
cussed in Sec. IV and summarized in Sec. V. In the 
Appendix we obtain upper bounds on some weighted av- 
erages of the initial amplitudes of the marked and un- 
marked states which are needed in the analysis. 



II. THE GENERALIZED GROVER ALGORITHM 

In the generalized algorithm the initialization step is 
modified as follows. Instead of initializing the register 
according to step 1, any initial distribution of marked 
and unmarked states can be used (e.g. the final state 
of any other quantum computation). In addition, the 
7r phase rotation of marked states described by IJ — 

^2x e i7rir ' x ) \x) (x\ is replaced by rotation by an arbitrary 
phase 7: 



(2.1) 



The rotation about the average, described by G = 
—WlfiW', where W is the Hadamard transform and 
I£ = I — 2 |0) (0|, is modified in two ways. First, the 
rotation of |0) by tt is replaced by the rotation of a prede- 
fined basis state \s) by an angle (3. Second, the Hadamard 
transform W is replaced by an arbitrary unitary opera- 
tor U. In this generalized algorithm, step 2B above is 
replaced by 



G 



where 



/f = /-(l 



)\s) 



s . 



(2.2) 



(2.3) 



In the generalized algorithm, the geometric interpreta- 
tion of step 2B as a rotation around the average ampli- 
tude of all states is not straightforward. We will demon- 
strate below that, in fact, one can return to this inter- 
pretation by i den tifying suit able variables. By inserting 
If from Eq. (O) into Eq. (Q we obtain that 



G = (1 - e*") \ V ) (v\ 



(2.4) 



where \rf) — U \ s). 



III. THE RECURSION EQUATIONS 



A. Analysis 



We will now analyze the time evolution of the am- 
plitudes in the generalized algorithm with a total of N 
states, r of which are marked. Let the marked amplitudes 
at time t be denoted by kiit), i £ M and the unmarked 
amplitudes by k (t) , i £ M where M is the set of marked 
states i\M\ — r) and M is the set of unmarked states 
(|M| = N - ry The initial amplitudes, fcj(O), i £ M 
and (0) , i £ M, at t = are arbitrary. Without loss of 
generality we assume that the number of marked states 
satisfies 1 < r < N/2. A general state of the system at 
time t will now take the form 



\9(t)) = J2 k i(t)\i) + ^2h(t)\i) 



(3.1) 



i£M 



A single Grover iteration Glj will transform the ampli- 
tudes kjit), j £ M to kj(t + 1) = (j\ Glj \g(t)) and the 
amplitudes lj(0), j £ M, to lj(t + l) = (j\ Glj \g(t)). We 
find that the recursion equations describing such itera- 
tion take the form 
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kj(t + 1) = (1 - e^ Vj + (1 - e^) Vj x 

ieM 

YjiiM - e^it) (3.2) 
ieTT 

lj{t + 1) = (1 - e^e^- £ ^(i)^ + (1 - e^jfe X 

5>(*H*-*i(*) ( 3 - 3 ) 



where 



(i \V) 



(3.4) 



If rji = for some i, the Grover iteration GJj changes 
only the phase of the state \i). Thus the probability to 
measure a state |i) with rji = remains constant. Hence 
we can treat the states with rji = separately from those 
with rji =/= 0. ^From now on we assume, without loss 
of generality, that for the given operator U, the prede- 
fined state | a) is chosen such that \rj) satisfies rji ^ 0, 
% = 1, . . . , N. We will now introduce new variables: 



Vj 



(3.5) 
(3.6) 



With these variables, the recursion equations will take 
the form 



k' 3 {t + 1) = (1 - e^K 7 E k 'M^ + C 1 " x 
ieM 

Yjmm?-^® (3.7) 

ieM 

l' 3 (t + 1) = (1 - e^e* X fcK*)M 2 + (1 " e i/3 ) x 

X^(*)l%| 2 -^(*)- ( 3 - 8 ) 



Let us define: 

2C(f) = (1 - e'V 7 X fc K*)l^i| 2 + (1 ~ ^) ]T 



ieM 



ieM 



(3.9) 



It becomes clear that the time evolution of all the am- 
plitudes (of both marked and unmarked states) can be 
expressed by 

kjit+l) = 2C{t)~e l ~*k' J (t), jeM (3.10) 
l' j {t + l)=2C(t)-l' j (t), jeM. (3.11) 

One can define the following weights 



ieM 



(3.12) 
(3.13) 



ieM 



quantifying the projections of the state \rj) on M and M, 
respectively. ^From the normalization of the state \rj) it 
is clear that W k +Wi — 1 . These weights can be used in 
order to define weighted averages of the new variables, of 
the form 



k'(t) = 
l'(t) = 



w k 

EieM^m 
Wi 



(3.14) 
(3.15) 



Using these averages to express C(t) we obtain 



2C(t) = (1 - e z/3 ) [e^W k k'{t) + W t l'(t)] . (3.16) 



By averaging over all the marked state s in E q. ( 3.10 ) and 
over all the unmarked states in Eq. ( |3 -11 ) we find that 
the weighted averages k'(t) and T'(t) obey the following 
recursion equations 



k / (t + l) = 2C{t)-e i ~<k'(t) 
l'(t + l) = 2C(t)-l'(t). 



(3.17) 
(3.18) 



These equations can be solved for k'(i) and l'(t), and 
along with the initial distribution this yields the exact 
solution for the dynamics of all amplitudes. We will pro- 
ceed to solve the recursion formulae for arbitrary complex 
initial amplitudes. Let us rewrite Eqs. (3.17) and (3. IS) 
in a matrix notation 



r(t + l) = Af(t) 



where 



and 



A 



k'(t) 
f(t) 



a b 



The matrix elements of A are given by: 



a = (1 
b=(l 
c=(l 
d=(l 



J0\ i-y 



y~<w k 

)Wi 



oil 



)e ll W k 



Wl - L 

The time-evolution of k'(t) and l'{t) is given by 
r(t) = A*f(0). 



(3.19) 



(3.20) 



(3.21) 



(3.22) 
(3.23) 
(3.24) 
(3.25) 



(3.26) 



In order to obtain explicit expressions for k'[t) and V[t) 
we consider the diagonal matrix: 
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A D = S~ X AS - 



A+ 
A_ 



(3.27) 



The eigenvalues A± of the matrix A are the solutions of 
det(j4 — XI) = 0. They can be expressed as 



where 



A± 



LO± = 7T H ± id 



(3.28) 



(3.29) 



and the angular frequency id is in the range < id < w, 
and satisfies: 

+ 1 — 1 

cosuj = Wkcos 1 —— - + Videos-— -. (3.30) 



The matrix S, which consists of the corresponding col- 
umn eigenvectors, takes the form 



1 



1 



A+— a X-—a 



(3.31) 



We will now apply Eq. ( 3.27 ) to reconstruct the matrix 
A, and compute A . A simpler expression for the time- 
evolution is obtained: 



f(t) = SAo'S^riQ), 



where 



A D * = 



Aj 







A_* 

The time dependence is given by 

k'(t) = zie 4 ^' - z 2 e iu - t 
T{t) = z 3 e M+t - z A e luJ - f 

(A_ -a)k'{0)-bi'(0) 

A_ - A+ 
(A + - a)fc'(O) - bl'(0) 
A_ - A+ 



where 



Z2 



and 



A + — a 
Zz = ; Zi 



b 

A_ — a 



Zi = 



-Z2- 



(3.32) 



(3.33) 



(3.34) 
(3.35) 



(3.36) 
(3.37) 

(3.38) 
(3.39) 



We have now completed the solution for the time depen- 
dence of the averages k'(t) and However, our aim 
is to obtain the time evolution of the individual variables 
fcj'(i) and li (t) and from them to extract the amplitudes 
h(t) and U(t) . Subtracting Eq. j3T7|) from Eq. ( ^Tp| ), 
and Eq. fl3.18| ) from Eq. (|3.1l[ ) one finds that 



ki(t + 1) - k'(t + 1) 
Jj(t + 1)-P(t + 1) 



- e *r[*5(t)-fc'(i)] (3.40) 
-K(t) -?(«)], (3-41) 



namely, the difference, in absolute value, between each of 
the variables fe/(t), h'(t) and the averages of the corre- 
sponding sets are time independent. This means that 



Ak> = k>(Q) - fc'(O) 
A^(0)-f(0) 



(3.42) 
(3.43) 



are constants of motion. Thus, the time dependence of 
the variables follows 



k[(t) = k'(t) + (-l)'e^'A^ 
l' i (t) = 7(t) + (-l) t Al' i . 



(3.44) 
(3.45) 



The time evolution of the amplitudes can now be ob- 
tained: 



k i (t)=r H [k'(t) + (-l) t e i i t Ak' i ] 
h(t) = '(*) + (-1)* Al'A. 



(3.46) 
(3.47) 



In this picture all the marked as well as the unmarked 
states evolve in unison so it is sufficient to follow the 
time evolution of the average in each set. The only fea- 
ture distinguishing the states from one another is their 
initial deviation from the average. 



B. Results 



^From Eqs. (3.44) and (3.45) it follows immediately 
that the weighted variances 



ieM 



(3.48) 
(3.49) 



ieM 



are time-independent and therefore at any time t they 
can be replaced by o\ = er|(0) and of = of(0), re- 
spectively. When a measurement is performed at time 
t, the probability that a marked state will be obtained 
is P(t) = J2ieM l^iWI 2 - Since all the operators used are 
unitary, the variables k/(t) and l/(t) satisfy the normal- 
ization condition 



Ei^ 2 i^)i 2 + Ei^i 2 i^)i 2 



i 



(3.50) 



ieM 



ieM 



at all times. Using the definitions of k^it) and an d 
t heir weigh ted a ve rages k'(t) and l'(t) given in Eqs. ( ]3~5| )- 
(|3~6l) and (|3.14|) - (|3.15|) , one can bring Eqs. (|^|) and 
(3.49|) to the form: 



E l**(*)l s 



W k a 2 k 



W k \k'{tt 



(3.51) 



ieM 
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El 

ieM 



h(tt 



Wta? 



(3.52) 



The first equation provides the probability to measure 
a marked state at time t, while the probability to mea- 
sure an unmarked state is given by the second equation. 
We will now try to examine the probability that a mea- 
surement at time t will yield a marked state, and its 
dependence on the rotation angles /3 and 7, the unitary 
transformation U, the predefined state \s) and the av- 
erages and standard deviations of the initial amplitude 
distributions of the marked and unmarked states. Using 



Eq. (3.34) with zi = \z x \e 1 ^ and z 2 = Ne^ 2 , the prob- 



ability that a measurement at time t will yield a marked 
state is given by a sinusoidal function of the form 



P(t) = P av ~ APcos[2(wt ■ 



(3.53) 



where 



AP = 2W k \ Zl \\z 2 \ (3.54) 
is the amplitude of the oscillations, tt/lu is their period, 
Pav = W k (\ Zl \ 2 + \z 2 \ 2 +a k 2 ) (3.55) 
is the average, or reference value of the probability, and 



20 = 01 - 4> 2 



(3.56) 



is the phase. These parameters are found to depend on 
the unitary operator U which is used in the algorithm, 
the predefined state \s) and the angle (5 by which its phase 
is rotated as well as the angle 7 by which the phases of 
the marked states are rotated. The dependence on the 
initial amplitudes enters only through the weighted aver- 
ages of the variables fe/(0) and 1/(0) and their standard 
deviations. 

It is thus observed that P(t) does not depend on any 
higher moments of the initial amplitude distribution. 
This is due to the fact that all the transformations in the 
quantum algorithm are linear. Therefore, ki(t) and U(t) 
can be expressed as some linear combinations of hi (0) and 
U(0). The only nonlinearity appears in the expression of 
P(t) as a sum of squares of the amplitudes of the marked 
states at time t. Therefore, powers higher than quadratic 
are excluded in the expression for P(t). Moreover, P(t) 
does not depend on any other linear combinations of the 
first and second powers of the initial amplitudes, except 
for the particular weighted averages that compose the 
first and second moments. This results from the fact 
that Grover's iterations maintain a large number of con- 
served quantities, particularly the variances <r? and a 2 
which are constants of motion. As a result, the time evo- 
lution of the amplitudes can be fully described by the 
time dependence of the averages k'(t) and l'(t). Due to 
the linearity of the transformations, k'(t) and l'(t) can be 
expresses as linear combinations of k'(0) and l'(0). The 
probability P(t) of measuring a marked state at time t is 



given by Eq. (3.51) in which the first term on the right 
hand side, which includes the second moment is a con- 
stant of motion. Therefore, no other quadratic forms can 
appear. The second term depends on the first moment 
at time t, which is related to the first moment at t=0 
through the recursion equations. Thus, the dependence 
of P(t) on the initial amplitudes is only through the first 
and second moments of their distribution. 



IV. DISCUSSION 

In order to examine the performance of the generalized 
Grover algorithm we will now evaluate the highest pos- 
sible probability P max = P av + AP that a measurement 
will yield a marked state. We will also evaluate the op- 
timal number of iterations T after which the probability 
Pmax is achieved. 

The limit of difficult search problems is obtained when 
N ^> r > 1. This is reflected in the fact that in the 
original Grover algorithm the probability of measuring a 
marked state immediately after the initialization step is 
Wk — r/N. The assumption that W k <§C 1 carries over 
to the generalized case discussed here. It is satisfied in 
all cases except for very unlikely choices of the state |s) 
from which one Grover iteration with the operator U is 
sufficient in order to measure a marked state with high 
probability. In the analysis below W k will be considered 
as a small parameter. The highest possible value of the 
probability to measure a marked state is: 



Prr 



W k (\ Zl \ + \z 2 \) 2 + W k a k 2 



(4.1) 



Consider the parameters z\ and z 2 that express the 
dependence of P max on the initial amplitudes and the 
phase rotation angles /3 and 7. The expressions for z\ 
and z 2 in Eqs. (3.36) and (|3.37| ) include (A_ — A + ) in the 
denominator. Using Eq. (|3.28|) it can be written as: 



A_ — A+ = 2ie l 2 sinw. 



(4.2) 



Expanding sin u> in powers of W k <C 1 , one finds that in 
case that the angles /3 and 7 are different, and the dif- 
ference between them satisfies \/3 — j\ — 0(1), (e.g. in 
radians): 



A_ - A 1 = 2^e ^ 2 sm if— 



71 



0(W k ), 



(4.3) 



namely, the denominator, in absolute value, is typically 
of order unity. In case f3 = 7 

A_ - A+ = 4ie lf3 sin^IU fc ^ + 0(W k ?) (4.4) 

and the denominator is of order y/Wk < 1 In the fol- 
lowing we discuss the search problem in these two cases 
separately. 
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A. Different Rotation Angles: j3 7^ 7 

Here we consider the case when f3 ^ 7 and the dif- 
ference between them is fixed and finite, and satisfies 
\/3 - 7I = 0(f). In this case, using Eqs. (|3~22|), (f3~23|), 



(1 



28) and the assumption that W k <C 1, the matrix ele- 
ments can be approximated by 



*7 



0{W k ) 



b= l-e if} + 0(W k ). 
The eigenvalues are 



J0 



o J 7 



0(W k ) 
0{W k ) 



when /3>7(0</3<27r, 0<7< 2ir), and 



A 1 = -e 1 



0(W k 



X- = -e lf3 + 0(W k ) 



(4.5) 
(4.6) 



(4.7) 
(4.8) 



(4.9) 
(4.10) 



when 7 > /3. In the Appendix it is shown that the initial 

— 1/2 

amplitude distribution satis fies \k (0)| = 0(W h ) and 

|P(0)| = 0(1). /From Eqs. §|g ) and (H|) we obtain 
that |«i I = 0(1) and \z 2 \ = 0(W k ~ 1/2 ), in case /3 > 7, 
while in case 7 > /?, |zi| = 0(W k ~ 1/2 ) and |z 2 | = 0(1). 



Therefore, in both cases, using Eq. (3.54) one obtains 



AP 



0{W}' 2 ) 



(4.11) 



namely, the amplitude of the oscillations is negligible. 
Thus, the probability to measure a marked state after 
any number of iterations cannot be signifi cantly larger 
than the probability, given by Eq. ( 3.51 ), to measure 
a marked state at time t = 0. We conclude that in 
this case the algorithm fails to enhance the probability 
of measuring a marked state. The angular frequency is 
u) = |/3 - 7I/2 + 0{W k ). Clearly, for such a high fre- 
quency, for which the period is of the order of only few 
steps, it is hard to exploit the oscillations since measure- 
ments can be taken only in discrete times and are likely 
to miss the highest point. The analysis presented above 
applies as long as W k -C (/3 — r y) 2 . The conclusions are 
in agreement with the results of Refs. p2],p5j. 



B. Identical Rotation Angles: (3 = 7 

In this case the matrix elements can be approximated 
according to 



a 
b 



1 - e 



0(W k ) 
0(W fe ) 



A± 

This gives rise to 



2e^sin^ 



(4.12) 
(4.13) 

W k i+0(W k ). (4.14) 



Zl,2 = -jWk 



iV(0)e l%b ±W fc 'jfc'(0) +0(1) (4.15) 



where i\> = (tt - /3)/2. Inserting Eq. (DJ) into Eq. (O]) 
we obtain 



P„ 



il'(Q)e l v + W** #(0)1 + |<l'(0)e** - W k *k'{0)\ 



W k a k 2 + 0(W fc ). 



(4.16) 



To simplify this expression we will now use the identity 
(|o + b\ + \a- b\) 2 = 2{\a\ 2 + \b\ 2 + \a 2 - b 2 \), where a 
and b are complex numbers. We will also replace Z'(0) by 
WjH'(O) [note that = 1 + 0{W k )\ and find that 



Pn 



1- w m 2 -- Wi|F(0)| 2 



^ fc |fc'(o)l 2 



1 

+ 2 



Wi|P(0)| 2 e 2iW+Q '- afe:) 
O(Wfc) 



l^|fc'(0)| 2 

(4.17) 



where jfc'(O) = |fc'(0)|e lQfc and i'(0) = U'£0)|e iQi . This 
is in agreement with the results of Ref. |20|], where the 
case P — 7 = tt was studied using a different approach. 
We observe that P max depends on the statistical prop- 
erties (averages and variances) of the initial amplitude 
distribution of the marked and unmarked states. For a 
given distribution, the probability of measuring a marked 
state is bounded by P max = \ — Wio 2 . This upper 
bound is reached when ip + a; — a k — 0, as well as when 
I' = or k' = 0. This optimization can be achieved by 
an adjustment of the the rotation phases to the value 
(3 = it — 2(a k — oti). However, this requires one to know 
the difference between the phases a k and a/, which is 
not generally available for an arbitrary initial amplitude 
distribution. 

The optimal case of P max = 1 can be obtained by using 
the predefined state \s) and applying on it the operator 
U, to generate the initial amplitude distribution. In this 
case the initial unmarked state variance is 07 = 0, the 
weighted initial averages are fc'(O) = 1 and l'(0) = 1, 
and the phases a k = ai = 0. Thus, executing the gen- 
eralized Grover iterations using the same unitary oper- 
ator U, the predefined state \s) as the initial state and 
a rotation phase of f3 = 7 = tt enables one to measure 
a marked state with the optimal probability P max = 1. 
The original Grover algorithm is a special case in which 
U is the Hadamard operator and the predefined state 
|,s) = |0...00). 

Consider an arbitrary initial distribution of r marked 
and N — r unmarked states, with known averages k(0) 
and Z(0) respectively. The probability P(t) that a mea- 
surement at time t will y ield a marked state is a sinusoidal 
function, given by Eq. ( |3~53| ). The highest value of P{t) 
is obtained at time T, for which the argument of the co- 
sine function satisfies 2{loT + 4>) = tt. Thus, the number 
of iterations T which gives rise to the highest probability 
of finding a marked state upon measurement is 
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T = 



ir-2<fi 
2uj ' 



where the angular frequency is 



uj — 2 sin — v W k 



0{WI' 2 ) 



(4.18) 



(4.19) 



An interesting case is the one in which the average 
and variance of the initial amplitude distribution are not 
known, but different runs of the algorithm use initial am- 
plitudes drawn from the same distribution. Naively, one 
could pick a random number of iterations T r and thus 
find a marked state with probability P(T r ). Correspond- 
ingly, the expected number of repetitions of the entire 
algorithm using the same T r would be l/P(T r ) until a 
marked state is found. However, P(T r ) could be very 
small . A b etter strategy is now shown. ^From Eqs. ( 3.3p| ) 
and (3.53) it follows that the period of oscillation of P(t) 
depends only on the unitary operator U , the predefined 
state \s) as well as on the rotation phase f3 = 7 used in 
the algorithm. The phase <j) depends on the initial ampli- 
tude distribution and is thus unknown. Consider the case 
where one runs the algorithm twice, taking measurements 
at times 2\ and T 2 respectively, where Ti — T% = n/ (2u>). 
^From Eq. ( |3.53 ) it is clear that in one of the two mea- 
surements the cosine expression will be negative so that 
P(T) > P av . Since the probability P(t) must be non- 
negative at anytime, P av > AP. Since P max = Pro+AP, 
we also find that P av > P max /2. Therefore, in one of 
the two measurements one obtains P(T) > P max /2. In 
this case one needs twice as many repetitions to obtain 
at least half the success probability compared to the case 
when the optimal measurement time is known. The slow- 
down is thus at most by a factor of 4. 



V. SUMMARY 



marked and unmarked states. The time T and the prob- 
ability P max also depend on the unitary operator U which 
is used in the algorithm, the predefined state \s), the an- 
gle (3 by which its phase is rotated as well as the angle 
7 by which the phases of the marked states are rotated. 
Moreover, it was found that in order for the algorithm 
to apply the two rotation angles must be equal, namely 
= 7- 
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APPENDIX A: 

In this Appendix we obtain upper bounds on |fc'(0)| 
and I £'(0)| in th e in itial a mpli tude distribution. Ac- 
cording to Eqs. (3.5), and ( 3.14 ) the initial distribution 
weighted average of the marked states is: 



fc'(0) = 



(Al) 



^From normalization it is clear that |fcj(0)| < 1 for any 
marked state. Therefore 

w k - w k - w k y ' 

where r is the number of marked states and \rj\ = 
max ieM \rji\. Since I77I 2 < J^ieM l^l 2 = w k one obtains 



In this paper we have generalized the recursion equa- 
tion analysis of Grover's quantum search algorithm pre- 
sented in Ref. We applied it to the large class of 
Grover's type algorithms in which the Hadamard trans- 
form is replaced by any unitary transformation and the 
phase inversion is replaced by a rotation by an arbitrary 
angle. We derived recursion equations for the time evo- 
lution of the amplitudes of the marked and unmarked 
states, for any initial complex amplitude distribution. 
These equations were solved exactly. ^From the solution 
we obtained an expression for the optimal number of iter- 
ations T after which the probability of finding a marked 
state upon measurement is the highest. The value of this 
probability, P ma x, was also obtained. Both T and P max 
are found to depend on the averages and variances of 
the initial amplitude distributions of the marked and un- 
marked states, but not on higher moments. This is due 
to the linearity of the transformations and the large num- 
ber of conserved quantities, particularly the (weighted) 
variances of the distributions of the amplitudes of the 
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W(0)\<^- = rW k 



(A3) 



Therefore, typical initial amplitude distributions in large 
search problems satisfy 



\k'(o)\ = o(w k -*). 



(A4) 



According to Eqs. ( p.q ) and ( 3.15 ) the weighted average 
of the initial distribution ^(0) satisfies 

p ( o) = E^'W) = J2 m*im + o{w k ). (A5) 



w, 



ieM 



Using Eq. (3.1) for the initial state of the system 
\g(0)}, the expression for fc'(0) [given in Eq. (Al)] and 
rji* = (s\ U* \i), one can bring Eq. (A5) to the form: 



P(0) = (s| U* \g(0)) - W k k'(0) + O(Wfc) 



(A6) 



Using Eq. (A4) one obtains 
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V(0) = {s\U* \g(0)) + O(W k -*). 

Since |<?(0)) is normalized and U* is unitary. 
(s| U* \g(0)) | < 1. Therefore, 



namely 



|P(o)| <i + o(w k - 
|F(0)|=O(i). 
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